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Motivated by recent experiments on Yb-doped CeCoIn 5 , we study the effect of correlated disorder in Kondo 
lattice. Correlations between the impurities are considered at the two-particle level. We use mean-field theory 
approximation for the Anderson lattice model to calculate how the emergence of coherence in the Kondo lattice 
is impacted by correlations between impurities. We show that the rate at which disorder suppresses coherence 
temperature depends on the length of impurity correlations. As impurity concentration increases, we generally 
find that the suppression of coherence temperature is significantly reduced. The results are discussed in the 
context of available experimental data. 
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I. INTRODUCTION 

Physical effects driven by an interplay of induced disor- 
der and strong correlations in condensed matter systems have 
been one of the central topics of experimental and theoret- 
ical research for many years. One of the active research 
directions is particularly focused on the properties of disor- 
dered heavy-fermion materials. In these materials, strong 
coupling between magnetic and non-magnetic degrees of 
freedom not only leads to a significant enhancement of the 
electron's mass below certain temperature T*, but also to a 
number of remarkable phenomena 1 , such as unconventional 
superconductivity 2 5 , magnetic-field-induced non-Fermi liq- 
uid metallic phases 6 , hidden-order phase transition 7 and, re- 
cently proposed topological insulating phases 8 . Studying how 
these systems react to disorder may help to better understand 
the microscopic mechanisms responsible for the formation of 
coherent states. 

Earlier experiments have been mostly focused on the ef- 
fect of doping on superconductivity and the onset of coher- 
ence (i.e. formation of heavy fermions) in Kondo lattices. 
It was generically found that a substitution of magnetic sites 
with non-magnetic impurities yields a substantial reduction in 
both coherence and superconducting critical temperatures^^]. 
Later similar studies have been performed on materials be- 
longing to the group of '115' heavy-fermion superconductors 
and have given qualitatively similar results 14 15 to what has 
been observed in other heavy fermions. 

The problem of how weak disorder affects the onset of co- 
herence in Kondo lattice has been discussed theoretically in 
Ref. Ifl3l . It was assumed that impurity atoms are indepen- 
dent and, therefore, the disorder correction to the conduction 
and /-electron self-energies were computed exactly within 
the self-consistent Born approximation. In particular, it was 
shown that impurities considerably affect the formation of the 
heavy fermion state when the strength of the disorder poten- 
tial is comparable to the single site Kondo temperature. The 
dominant contribution to the suppression of T* comes from 
the disorder potential on /-sites. It was pointed out in Ref. 
lfT3l . this is not surprising given the fact that the heavy quasi- 
particle spectral weight is dominated by the /-states. Overall, 
the results of this theory 13 were found to be in agreement with 
available experimental data. 



Recent measurements^^, however, consistently demon- 
strated that the substitution of Yb on Ce sites in CeCoIns 
does not lead to significant changes in the properties of this 
material. Specifically, according to the data reported in Ref. 
lfT6ll . the onset of coherence remains essentially unaffected 
with Yb-doping, while the critical temperature of the super- 
conducting transition decreases linearly with Yb concentra- 
tion. Given the presence of strong correlations in this mate- 
rial, manifested by the linear temperature dependence of re- 
sistivity at low temperatures 18 and unconventional af-wave su- 
perconducting pairing, these experimental results are highly 
unusual. 

The recent experimental data motivate us to revisit the prob- 
lem of disorder in heavy fermion materials. Specifically, in 
this paper we will address the following question: how the 
onset of coherence in the Kondo lattice as a function of impu- 
rity concentration is modified by the correlations between the 
impurities? To answer this question, we will adopt the theory 
developed in Ref. Ifl3l appropriately modified to include the 
correlations between impurities. Specifically, impurity corre- 
lations are described by the probability distribution function 
which depends on the distance between the two impurities. 
Within our theory we will show that even in the presence of 
weak disorder (i.e. when the self-energy corrections are small 
compared to the relevant energy scales in the system), the ex- 
istence of the pair impurity correlations significantly slows the 
suppression of the coherence temperature. We also draw mo- 
tivation from the experimentally observed enhancement of the 
superconducting critical temperature in high-7\. superconduc- 
tors caused by the ordering of the oxygen interstitial*^, as 
well as from a weak suppression of the superfluidity in 3 He in 
the aerogeP^I 

The correlations between the impurities in a material can, 
in principle, be induced by the local lattice strains. To un- 
derstand why the impurity correlations may slow down the 
suppression of the coherence temperature, let us first intro- 
duce the characteristic length scale R on which impurity dis- 
tribution function significantly deviates from unity. Then, 
within the Born approximation, one can show that there will 
be two contributions to self-energy. One contribution, Eg, cor- 
responds to the scattering of electrons on the same impurity 
and, upon the averaging over disorder, this contribution is pro- 
portional to the concentration of impurities n,- mp . The second 
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contribution, Ey, describes the scattering of electrons on two 
different impurities and, therefore, is proportional to nf mp . In 
the presence of impurity correlations, however, Ey becomes 
proportional to nj mp R 3 . Thus if the radius of correlations is 
large enough (i.e. rij mp R 3 ~ 1), becomes comparable with 
the first, linear in n lmp , self-energy correction Whether Sy 
can significantly compensate for Z„ will, of course, depend on 
the specific form of the impurity distribution function, p(r/R). 
For the specific form of p(r/R), we will show that increase in 
R yields higher values of T* compared to ones obtained for 
the case of uncorrelated disorder. 

Our paper is organized as follows. In Section II we dis- 
cuss the possible origin of the impurity correlations in heavy- 
fermions and the self-energy corrections resulting from corre- 
lated disorder. In Section III we present the calculation of the 
coherence temperature with the slave boson mean-field the- 
ory in the presence of correlated disorder. In Section IV we 
provide a general discussion of our findings in the context of 
available experimental data. Our results are summarized in 
Section V. Finally, in Appendix A we discuss how the spe- 
cific choice of the impurity correlation function influences the 
electron scattering lifetime. 



II. PRELIMINARIES 



Generally speaking, the local ionic vibration of the impurity 
atom causes the change in A,. Expanding the energy splitting 
parameters up to the linear order in lattice displacements u,, 
we find 
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where A,o is the value of the energy splitting at equilibrium, 



is a strain tensor at site i 
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and the components of y y i define the deformation potential. 
Consequently, the lattice deformation at site i may effect the 
impurity at site j. This can be shown by formally integrating 
out the fields k® in the path-integral formulation of the prob- 
lem. Alternatively one can use the condition for the balance 
between the internal and external stress in the system, which 
yields^ 



H imp = 2^ A ;0 o-,v + 2_j u< 



(4) 



In this Section we define the impurity distribution form fac- 
tor: the central quantity for our subsequent analysis of the co- 
herence in the disordered Kondo lattice. In this Section for 
simplicity we limit the discussion to the problem of conduc- 
tion electrons moving in the disorder potential, which is as- 
sumed to be weak. 

We will first demonstrate that the interactions between the 
impurities may be induced by the local strains in the host ma- 
terial. Then we proceed with evaluation of conduction elec- 
tron self-energy, which will involve the impurity distribution 
form factor. Impurity distribution form factor formally ap- 
pears due to averaging over the position of impurities. The 
form factor is, in turn, determined by the impurity correlation 
function. We expect the impurity correlation function to be 
a functional of the interaction potential between impurities, 
which allows us to approximately determine the shape of the 
former. 



A. origin of the impurity correlations 



The interaction potential decays as C7y ~ r.r at large dis- 
tances. Most importantly, the sign of the interaction potential 
is defined by the local strain fields and can be either negative 
or positive, which corresponds to effective attraction or repul- 
sion between impurities. 



B. disorder distribution: form factor 

To make our discussion self-contained, in this subsection 
we briefly review the theory of correlated disorder first dis- 
cussed in Ref. 11221 for the problem of the superfluid 3 He in the 
porous medium. Hamiltonian describing the electrons moving 
in the disordered lattice reads: 



« £ J" ^ 3 r^<J(r-R,)^(r)^(r), 



(5) 



The two-impurity correlation can in principle be induced 
by the local lattice strains. To show this, let us first assume 
that the single impurity can be in only two states. The energy 
separating the these two states on impurity at site i is A,. Then 
the impurities are described by the following Hamiltonian: 



H, 



imp 



(l) 



where cr iz is a third Pauli matrix and summation goes over the 
impurity sites. 



where u is a strength of the disorder potential, (Ai-(r) creation 
operator of an electron with spin projection cr at point r in 
space, £(V) = -V 2 /2m - fi and fi is a chemical potential. 
For weak disorder, the second order correction to the elec- 
tron's Green function @(v, r'; r) = -(J T \/i lT (v, r)^o-(r' 0)) in 
momentum and Matsubara frequency representation is^ 



£ (2, (k,k';/w) = |wr£r w (k; ico)& m (k'; iai) 
X V f i!^ e -'(k-q)R.-Xq-k')R^(0) (q;/w) 



(6) 
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To recover the spatially homogeneous expressions fo r the 
correlation functions, we follow the standard procedur e 1 22 * 25 ! 
of averaging over the impurity positions R, : 



Y j(e i<*-vm )dis = nimp s(k-k'), 

i 

^-/(k-^RHtq-lOR,^ = n . mp6(k _ k0[1 + S(k _ q)]> 



(7) 



where we have introduced the impurity distribution form fac- 
tor S (k): 



S(k-q) 



L unp 



I 



tTrp(r)e 



/(k-q)r 



(8) 



Here n imp - Ni mp /V is the concentration of impurities and 
p(r) is the impurity distribution function, which depends on 
the distance between two impurities. Thus, for conduction 
electron self-energy we obtain: 



E(k, iof) — rij, 



[l+S(k-q)]g ((, )(q;i W ) (9) 



If we set the origin of the reference frame at one impurity, 
p(r) gives the probability to find a second impurity at distance 
r from the first one. At large distances p(r — > oo) — > p x . 
After the proper choice of the normalization constant we can 
set poc, — 1. Thus, for uncorrected imputiries p(r) = 1 and 
it follows 5(k) = rii mp 6(k). In particular, this implies that 
the second term in the second equation |7| is of the order of 
0(nf mp ) and therefore can be neglected. However, as we will 
show below, for the correlated disorder, characterized by the 
correlation radius R, it is possible that the form factor becomes 
of the order (9(1) and thus it needs to be taken into account. 

Since at large distances impurities are not correlated, ex- 
pression for the form factor S(k) contains both correlated 
and uncorrelated contributions. To single out the contribution 
from correlations we add and subtract the uncorrelated part: 



— ft imp ^ 

— ft imp ^ 



dMp(r) - l]e ,kr + n imp S(k) 



d 3 rv(r)e' kr + n imp 5(k), 



(10) 



where function v(r) serves as the measure of impurity corre- 



lations. We can also omit the last term in ( 10 1 since it yields 



a trivial correction of the order of 0(nf mp ) to the self-energy. 
Therefore, to evaluate the momentum dependence of the form 
factor, we first will have to specify the function p(r) or, alter- 
natively, correlation function v(r). 

From general considerations^it follows that the correlation 
function v(r) should be a functional of the impurity interaction 
potential I/y = t/(|r, -r ; |), Eq. Q. If we were to consider 
a weakly interacting ideal gas of atoms, for the pair of two 
atoms one could write the correlation function as Vfi(r, ; ) ~ 
[exp(-Uij/T) - 1], where T is a temperature (see Ref. l26i l. 

We will c onsider two correlation functions of the following 

forrrpmi 



y±(rij) = ± 



A ± 



R 



3-a 



- 1 



-ru/R 



(11) 



where 1 < a < 3, which also shows qualitatively similar de- 
pendence on distance as Vfi(r,j) (with corresponding attractive 
or repulsive potential U/j). Our choice of the correlation func- 
takes into account the possible fractal structure for 
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tions 

the impurity atoms by means of the fractional exponent a and, 
as it turns out, makes the effects of the impurity correlations 
more pronounced compared to the ones described by v B (rij). 
We allow for the possibility of having two types of correlation 
functions (Hi, which qualitatively correspond to an effective 



attraction or repulsion between the impurity atoms on inter- 
mediate distances. This is done in order to understand better 
the origin of the "healing" effect in Kondo lattices, which we 
will discuss in the next Section. 

To evaluate the value of the constants A+ in Eq. 



1 1 1, let 



us introduce the impurity density distribution function n(r). 
Consequently, the following correlation function is considered 
(for details, see discussion in Ch. XII of Ref. l26ll ): 

<«(!-! )n(r 2 )> = n imp 6(ri - r 2 ) + "Lp/^ri - r 2 ). 

For the fluctuations in the impurity density one finds 



{6n(r x )5n(r 2 )} = n imp 5(v\ - r 2 ) + nf mp v ± (r) 



(12) 



Integrating both parts in Eq. ( 12 1 with respect to the spacial 
coordinates and taking into account that in our case that the 
fluctuations in number of impurities are absent, we fincP 



4nn, 
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„„,, i v ± (r)r dr = -\. 
b 



(13) 



The integral ( 13 i with v ± (r) given by (Hi can be evaluated 



exactly. We find 



F(a) \ 4nn imp R 3 



The momentum dependence of the form factor ( 1 0| 1 1 1 can be 
evaluated exactly. We find 



S ± (k) = +Ann imp R 



3 , A ± F(a - 1) sin[(a - 1) tan" 1 ^)] 
kR[k 2 R 2 + 1](«-D/ 2 



[1 +(/W?) 2 ] 2 



(14) 



Note that the "strength" of the impurity correlations is effec- 
tively measured by the value of the dimensionless parameter 
ftimpR^ ■ From Eq. ( 14 1 we immediately observe that when the 



impurity correlation length scale R is large enough, the form 
factor becomes of the order 0(1). In Appendix A we evaluate 
the effect of the impurity correlations on the electron's scat- 
tering time, which is determined by the imaginary part of the 
self-energy. 

In the next Section we study how the correlated disorder 
will affect the onset of the coherence on the Kondo lattice: 
a periodic lattice of predominantly localized /-electrons hy- 
bridized with conduction electrons. 
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in. DISORDERED KONDO LATTICE 

To describe the physics of an interplay between the disorder 
and coherence in the Kondo lattice we will employ the Ander- 
son lattice model (ALM) in the limit of infinitely large Hub- 
bard interaction U between the localized /-electrons. Since 
the disorder is assumed to be weak, we first briefly discuss the 
slave-boson mean-field approximation for the ALM and intro- 
duce the coherence temperature associated with the formation 
of the heavy fermion metallic state. We than proceed with the 
discussion of the disorder effects on the value of the coherence 
temperature for both correlated and uncorrected disordei^. 



A. mean field theory of the Anderson lattice model 

We begin with writing down the model Hamiltonian to de- 
scribe the physics of the disordered Kondo lattice. The Hamil- 
tonian describing conduction electrons is 



(15) 



k (7=U 



where ^ = k 2 /2m - fi is the dispersion of conduction elec- 
trons (assumed to be parabolic), fi is a chemical potential, cr is 
a spin and c. is a conduction electron creation operator. Con- 
sequently, the Hamiltonian which describes the /-electrons is: 

Nr 

", ~ X e 'H J - + U Yj&f™flfi<s- (16) 

j,a= 1 iaa' 

where fl creates an /"-electron on site / in a state of a low- 
est lying multiplet AV-degenerate multiplet of the /-ion, et is 
the /-electron energy and U > is the strength of the Hub- 
bard interaction between the /-electrons. We note that index a 
is not a spin index due to the presence of the strong spin-orbit 
coupling. Generally states belonging to the multiplet T are de- 
scribed by the total angular momentum J and z-component M 
or some linear superposition of those states and in the second 
term ( fT6} the summation is restricted to a + a'. 

Finally the term describing how conduction electrons are 
hybridized with localized /-electrons is 



ZVr 



+ V*f] a c u 



(17) 



Here V is a hybridization matrix element between the con- 
duction electrons and localized /-electrons and ct are project 
the conduction states with spin cr onto an /-state classified by 
the z-component of angular momentum M. This projection 
is furnished by the presence of strong spin-orbit coupling 27 
and, generally speaking, also yields the hybridization matrix 
element to be non-local. However, as we have verified, the re- 
alistic spatial (or momentum) dependence of the hybridization 
elements does not affect our main results and, therefore, will 
be ignored here. Thus, in the absence of disorder the periodic 
Anderson model Hamiltonian reads: 



Hp 



= H C +H f + H, 



(18) 



In heavy-fermion materials Hubbard interaction is much 
larger than the bandwidth, so that we can take it formally to 
infinity, U — > oo. In this case, the multiply occupied states, 
f 2 , / 3 etc., on which Hpam ( 18 I operates must be projected 



ouPh This is achieved by introducing the slave boson field 



bi, which act like projectors^ES. As a result for our model 
Hamiltonian we can write: 



H PAM -II, - X [ Vc ]* b ]fj« + V ' n„<-i:l>j\ 



j,a=l 
j,a=l 



(19) 



It has to supplemented by the constraint condition of not al- 
lowing for more than one /-electron per site: 



A'r 



(20) 



To find the ground state governed by the model ( 19|20 i the 
mean-field approximation is adopted. Within the mean-field 
theory, the slave-boson operators are replaced with their aver- 
age values: 



h -> (bt) = a, 



(21) 



which is then needs to be determined self-consistently. The 
self-consistency equations are derived from minimizing the 
free energy with respect to the slave boson amplitude a and 



also Lagrange multiplier A used to enforce the constraint ( 20 1. 
In addition, the total number of particles needs to be con- 
served. As a result one obtains the following system of three 
nonlinear equation ^ 31 ! 34 ! 

(A-e f )a + VTY j fc (k,a> n ) = O, 

k,£j„ 

(« 2 -?r) + r^%(k ; w„) = 0, 
(q r - a 2 ) + T X @cc(K w„) = C, 



(22) 



where T is a temperature, qr — 1 /Nr and C is a constant (for 
heavy fermion metals C - 1/2). Lastly, the propagators en- 



tering into Eqs. (22 1 are given by: 



@cc(K a>„) = 
0/ c (k, w„) = 

@ff(K (O n ) = 



On ~ £f 



(io) n - e f )(i(o„ - £ k ) - V 2 a 2 ' 
Va 

(ia>„ - e f )(io) n - f k ) - V 2 a 2 ' 

ioJ n ~ & 

(ia) n - e f ){iu„ - £ k ) - V 2 ^ ' 



(23) 



The mean-field approximation becomes exact in the limit 
A^r — > oo while for the finite values of Nr the corrections 
to mean-field results are of the order of 0(1 /Nr). We note, 
however, that the slave-boson mean-field theory proved to 
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be very reason able approximation for many heavy-fermion 
material 9^21311 

The nonzero solution for the amplitude a appear for the first 
time at temperature 



Tq =a D exp 



k/l 



PfV 2 



(24) 



where D is a width of the conduction band and pf is a density 
of states at the Fermi level. Non-zero value of a signals the 
opening of the hybridization gap in the single particle spec- 
trum. Opening of the gap also yields the enhancement of 
the electron's masJ^l. On the other hand, experimental data 
shows that in heavy fermion metals resistivity has a maxi- 
mum at some temperature, which is usually associated with 
Tl, and then decreases upon further cooling 1 4 35 . Therefore, 
Tl is interpreted as a coherence temperature below which the 
heavy-fermion metallic state is formed. 



B. effect of disorder on coherence in Kondo lattice 

In this Section we discuss the effect of the correlated disor- 
der on the coherence temperature T*. Although we are well 
aware of the fact, that slave-boson mean-field theory may not 
be entirely reliable approximation to compute T*, we believe 
that it suffice to demonstrate disorder correlations affect the 
value of T* in comparison with uncorrected case. 

Since disorder breaks translational symmetry of the lattice, 

as well as other mean-field 



the slave-boson amplitude (21 



parameters become inhomogeneous. However, for the case 
of weak disorder, i.e. when the self-energy corrections to 
the conduction and /-electron propagators are assumed to be 
much smaller than T*, we can perform disorder averaging. 
This would yield the parameters of the theory homogeneous, 
so that to determine the value of T* in the presence of disorder 
we still need to solve the system of equations p2) with propa- 
gators including the self-energy corrections due to disorder 13 . 
The disorder is simulated by the following Hamiltonian: 



WmjfJj>jCja + h-C.)] 



(25) 

We assume that random variables W/,-, W c i and W„„- are 
Gaussian random variables, satisfying {Wf yC jn;i)dis = and 



Wi da, where (...),«. denotes the aver- 



aging over the random distribution of Wf >c< ^¥^. 

Disorder corrections to self-energy are most compactly 
written using the matrix notation for the electron propagators: 



QffQtJo) @ fc (k,iu)) 
c f(k,ia>) & cc (k, iu) 



(26) 



Within the Born approximation for the self-energy we find: 



/»3 
— ^W^ (q, MW [1 + S + (k - q)] . 
(iny 



where the disorder averaging has already been carried out, im- 
purity correlation function S +(k) is given by ( fT~4| > with a = 3/2 
and the elements of the matrix W are self-evident. Thus, the 



propagators entering into Eqs. ( |22| ) should be replaced with 

(28) 



G 1 - Qq 1 - E. 

In principle, we could have used the self-consistent Born 
approximation (SCBA) which amount to replacing the bare 



propagator in ( 27 1 with an exact one ( 28 i and solving for 
the self-energy self-consistently. This task can be easily ac- 
complished for the uncorrelated disorder when the self-energy 
does not depend on momentum. In the presence of the disor- 
der correlations, however, we have chosen the ladder approx- 
imation for the self-energy in order to facilitate the numerical 
calculations. We have also verified, that the self-energy com- 
puted within SCBA does not differ significantly from the ones 
given by p7| ) for a wide range of the disorder strength. In ad- 
dition, to keep our results tractable we will adopt the same ap- 
proximation as in Ref. lfT3ll and ignore the disorder potential 
in hybridization, W m = 0. We have checked, that this approx- 
imation does not lead to any qualitative changes in the values 
of a and T* for both correlated and uncorrelated disorder. 

We evaluate the coherence temperature T* for fixed dis- 
order potential W = Wf = W c = 0.85T*. This specific 
choice for disorder strength is dictated by the fact that both 
slave-boson amplitudes and coherence temperature should be 
significantly suppressed when disorder becomes of the order 
of the single ion Kondo temperature 13 . This, in turn, should 
make an effect of disorder correlations more pronounced. On 
Fig. [T]we present our results of the numerical solution of Eqs. 
d22b for T* /Tl as a function of disorder w — n lm „ W 2 for vari- 
ous values of the impurity correlation radius R/% (£ = Vf/TT). 




(27) 



FIG. 1: Dependence of the coherence temperature (in the units of 
the coherence temperature for the clean system) on impurity con- 
centration n imp multiplied by the square of the disorder potential, 
w = n imp W 2 . Normalization parameter w c corresponds to the crit- 
ical value of disorder when the coherence temperature vanishes in 
the absence of impurity correlations. 

From Fig. [T]we see that impurity correlations have small ef- 
fect at small concentrations, which is expected given the fact 
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the self-energy corrections due to scattering off two impuri- 
ties is proportional to « lmp /? 3 « 1. As H, m/ >/? 3 increases, the 
suppression of the coherence temperature becomes less pro- 
nounced. It is also important to remember that as the value of 
the parameter ni„, p R 3 reaches a certain critical value, the lad- 
der approximation, used to evaluate the correlation functions, 
breaks down. This is the reason why in Fig. 1 we have re- 
stricted the range of disorder parameter variations to the value 
w/w c = 1.25. We also find that the the same calculation with 
the correlation function v_(r), Eq. ( fTT) , yields even stronger 
suppression of T* as correlation radius increases, so that no 
"healing" effect occurs. 



IV. DISCUSSION 

The rate at which the coherence temperature is suppressed 
as the correlation radius is changing clearly depends on the 
specifics of the spacial dependence of the impurity correlation 
function v ± (r), Eq. (Ill 



The correlation function v + (r) at 
distances r > R is negative which can be interpreted as an ef- 
fective attraction between the impurity atoms. Thus we lead to 
interpret the observed enhancement (in comparison with un- 
correlated case) of T* as being due to the cluster formation of 
Yb atoms. In this sense our interpretation is therefore simi- 
lar to the one discussed in Ref. [17]. The difference, however, 
is that the electronic state, although being spatially inhomoge- 
neous due to disorder, is not strongly coupled to the impurities 
and the reduction in the T* suppression rate is a property of 
an impurity system itself: there is no special re-organization 
of the electronic system in response to disorder. 

In addition to insensitivity of T* to Yb-doping in CeCoIn 5 , 
equally puzzling property of this material is an unexpectedly 
slow suppression of the superconducting transition tempera- 
ture rj 16 ! 17 ! One intriguing question in this context is whether 
concentration dependences of both T* and T c are correlated 
with each other. We believe that this question can certainly 
be addressed within the recently proposed theory of the tan- 
dem pairing in CeCoInJ^D. In this theory, the Cooper pairing 
emerges as a result of strong hybridization between singly oc- 
cupied /-states and conduction electrons in the presence of 
two conduction channels emerging from the fluctuations into 
empty and doubly occupied /-states. The symmetry of the 
conduction channels is determined by the point group irre- 
ducible representations and, as a result, ultimately determines 
the symmetry of the superconducting order parameter. Within 
the large-N mean field theory, however, the emergence of the 
superconductivity corresponds to the opening of the second 
hybridization gap at the Fermi level. Therefore, within this 
picture, the behavior of both T* and T c will have to be corre- 
lated. We are planning to address this question in the future. 

Another important question is how the presence of correla- 
tions between the impurity atoms can be probed experimen- 
tally. We think that one way to probe impurity correlations 
is to use the scanning tunneling microscopy (STM) measure- 
ments. Using this technique, one should be able to measure 
the position of impurities and any non-trivial physics associ- 
ated with them. Alternatively, one can use momentum space 



probes such as X-ray diffraction spectroscopy-^ to resolve for 
the presence of any fractal structures in the array of Yb atoms. 



V. CONCLUSIONS 

In this paper we have revisited a problem of how disor- 
der affects the formation of heavy-fermions in Kondo lattice. 
Specifically, we have studied how the presence of correlations 
between the impurity atoms influences the value of the co- 
herence temperature. Our calculations have been performed 
within the framework of the slave-boson mean-field theory. 
Although this theory may not provide an entirely reliable esti- 
mate for the coherence temperature, it still allows us to study 
the effects of impurity correlations. In addition, we have also 
restricted ourselves to the ladder approximation to evaluate 
the disorder corrections to self-energy. This approximation 
limits the range of applicability of our theory to the case of 
strongly correlated impurities, but still is reliable when the 
correlations are not very strong. In agreement with previous 
studies we found that in the absence of correlations the coher- 
ence temperature is strongly suppressed with increase in im- 
purity concentration or values of disorder potential. When the 
correlation radius becomes comparable to the heavy-fermion 
coherence length we find the suppression of the coherence 
temperature is significantly reduced. Our theory can be di- 
rectly applied to the case of Yb-doped CeCoIns to account 
for the dependence of coherence temperature on Yb concen- 
tration. In particular, we suggest that the Yb impurities are 
strongly correlated in this material, which leads to insensitiv- 
ity of coherence temperature to disorder. The specific type of 
correlation potential can, in principle, be probed by the STM 
measurements. 
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Appendix A: electron scattering time due to correlated 
impurities 

In this Section we evaluate the changes in the time between 
two impurity scattering events due to the correlated disorder. 
We define the electron's scattering time due to disorder as fol- 
lows 



Im[E ± (* F ,0)]. 



(Al) 



where p F is a Fermi momentum. We evaluate t ± numerically 
using ([9]) and impurity correlation functions (111 at 7" — > 
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and fixed impurity concentration n„„ p = 1% (per unit vol- 
ume) as a function of the dimensionless parameter nj mp R 3 . 
We show the results calculation on Fig [2] As we can see 
from Fig. [2] depending on the sign of the correlation function 
v ± (r), electron's scattering time is either drastically increased 
or discreased with an increase in impurity correlations until 
impurity correlation parameter g ~ 1 . For higher degree of 
correlations, the changes in scattering time are moderate. 



FIG. 2: Plot of the electron's relaxation times t ± , Eq. | |A1} , for cor- 
related weak disorder as a function of the dimensionless parameter 
n jmp R 3 for two values of the exponent a, Eq. (TTJl. We use the col- 
lision time To for uncorrected disorder to set the scale. As we have 
expected, the correlations between the impurities lead to an increase 
between the electron impurity scattering events. 




Our results in this Section demonstrate how the presence of 
correlations between impurities affect the scattering properties 
of conduction electrons. By comparing the signs of the func- 
tions Vfi(r, 7 ) and v ± (r (/ ), Eq. ( [TTj i, we can associate the sign of 
the impurity interaction potentials with the sign of the corre- 
lation function. In particular, it follows that the interaction de- 
scribed by the correlation functions v+(r) are repulsive at short 
distances r «c R and attractive at large distances r » R, in- 
teractions v_(r) are repulsive for all distances. Our results can 
be qualitatively understood as follows: the repulsion (attrac- 
tion) between impurities means a larger (smaller) mean-free 
path for conduction electrons and, as a consequence, increase 
(decrease) in scattering time. Independent of the sign of the 
form factor, however, the effect of impurity correlations is the 
renormalization of the total scattering cross section. 
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